Excitons in Graphene and the Influence of the Dielectric Environment 
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The exciton Wannier equation for graphene is solved for different background dielectric constants. It is shown 
that freestanding graphene features strong Coulomb effects with a very large exciton binding energy exceeding 
3 eV. A second-order transition to a weak Coulomb regime is found if the effective background dielectric con- 
stant exceeds a critical value. All bound-state solutions vanish for epitaxial graphene on a substrate with large 
background permittivity, such as SiC. 



Graphene is an effectively two-dimensional (2D) system 
with unusual electronic and optical characteristics. It is widely 
believed that the electronic properties result from the sp 2 hy- 
bridization between one s and two p orbitals of the constitut- 
ing carbon atoms, leading to the typical planar honeycomb lat- 
tice. The remaining p orbital forms a half-filled it band respon- 
sible for the electronic and optical properties ilHj. Using the 
tight binding(TB) approximation with nearest neighbor hop- 
ping only, it was first shown by Wallace |1[| that the hexagonal 
symmetry of the honeycomb lattice leads to the formation of 
symmetric electron and hole bands that touch each other at 
two nodes. In the vicinity of these so-called Dirac points, the 
dispersion is linear and the density of states vanishes. Accord- 
ingto this picture, which is supported by many experiments 
|4|-|9j], graphene behaves as a semimetal whose massless Dirac 
electrons have a Fermi velocity Vp ps 10 6 m/s. 

Whereas the existence of quasi-relativistic electrons is gen- 
erally accepted, much less agreement exists concerning the 
role and the con seq uences of Coulomb interaction effects 
in graphene I10h25I1 . Studies based on a renormalization- 
group analysis predict a stable semimetallic ground state 
1 10H13I1 . however, nonperturbative methods yield a semimetal 
to excitonic insulator transition at sufficiently high coupling 
strengths lfl4l-[l8l. 



As a convenient measure of the relative strength of the 
Coulomb interaction, one can use the effective fine-structure 
constant «c = e 2 /4ireoehvF, where e = e r or (1 + e r )/2 is 
the effective background dielectric constant for a single layer 
of graphene embedded in, or grown on top of a substrate, re- 
spectively. As we discuss in this Letter, the Coulomb inter- 
action introduces a sensitivity to the dielectric environment 
via background screening, contradicting the widely believed 
paradigm that the physical properties of graphene are largely 
independent of the environment. For freestanding graphene in 
vacuum «g ~ 2.41, indicating prominent Coulomb interac- 
tion effects. 

It is well known from semiconductor physics that the 
Coulomb attraction between electrons and holes may lead 
to the formation of bound electron-hole pairs, i.e. excitons. 
While excitonic effects in metals are generally believed to be 
of minor importance because of strong screening, in graphene 
the screening of the long-ranging part of the Coulomb inter- 



action is suppressed by the vanishing density of states at the 
Dirac points. Hence, it is not surprising that excitonic reso- 
nances have been observed in one-dimensional metallic car- 
bon nanotubes M26I1 . Excitonic binding energies in these sys- 
tems ranging from 50 to 100 meV have been calculated by 
solving the Bethe-Salpeter equation and within a density ma- 
trix approach 122112711 . For planar 2D graphene, first-principles 
calculations have become available only recently 120(1 . They 
predict a self-energy correction of the Fermi velocity that is in 
good agreement with experimental findings |5_I| and an exci- 
tonic shift of the dominant optical absorption peak as large as 
600 meV. Hartree-Fock based configuration-interaction calcu- 
lations for graphene quantum dots yield an excitonic redshift 
of approximately 300 meV 11911 . 

In this Letter, we show that the relative importance of the 
Coulomb interaction effects in graphene is dominated by its 
dielectric environment. In particular, we predict that free- 
standing graphene in vacuum has Coulombic properties that 
are very different from those of epitaxially grown graphene 
on substrates like Si02 or SiC with relative e in the THz range 
of 4.45 and 9, respectively. 

A useful and convenient criterion for the importance 
of Coulomb interaction effects is the existence of bound 
electron-hole pairs, i.e. excitions. The prediction of a finite 
exciton binding energy in a gapless system like graphene, is 
an indication for the fact that the system energy an be reduced 
below the tight binding (TB) ground-state level by exciton for- 
mation. In order to investigate this feature, we consider the 
infinitesimal transformation U([i) = exp(/3B^) generated by 
the exciton operator B* x = ^ k ^x(k)el/il k . where <f>\(k) 
is a normalized wave function, and e k and h* are electron 
and hole creation operators defined within the bands. Start- 
ing from the TB ground state, an infinitesimal transformation 
yields an energy shift 
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Here, £k = £k + £ k are tne electron and hole dispersion re- 




FIG. 1. Coulomb matrix element for one sheet of graphene. The 
full Coulomb potential (solid line) and the Ohno potential (short- 
dashed gray line) are shown together with two asymptotic forms 
(long-dashed and dotted lines). Inset: Band structure and a single 
electron-hole-pair excitation in graphene. 



lations within the TB approximation yj-yl] and \\ are the 
Coulomb matrix elements calculated with the electron and 
hole wave functions 11281 12911 . 

The shift SE is negative if the graphene Wannier equation 
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has bound state solutions with E\ < 0. Eq.|2]is a direct gener- 
alization of the well known Wannier-equation in semiconduc- 
tor physics, and describes the relative motion of a Coulomb- 
interacting electron-hole pair with zero center of mass mo- 
mentum (c.o.m)[30]. Altough a general dispersion does not 
allow for the decoupling of the relative and center-of-mass 
(c.o.m.) motion problems 13 ll 13211 . the smallest pair energy is 
obtained for states with vanishing c.o.m. momentum, as the 
kinetic c.o.m. energy is always positive. Hence, bound state 
solutions of Eq. |2] clearly limit the validity range of a weakly 
interacting model system. Within the linear approximation, 
the graphene Wannier equation is equivalent to the Dirac two- 
body problem with zero c.o.m. Recent real space studies con- 
sidering the Dirac Coulomb problem predict instabilities if the 
effective fine structure constant exceeds a critical value 1132 - 
3511 . Although restricted to positiv energies, the adressed top- 
ics are closely related to the exciton problem and indicate a 
nontrivial dependence on the background dielectric constant. 
To characterize the influence of the dielectric environment, we 
calculate the binding energy of electron-hole pairs with zero 
c.o.m. as a function of the effective background e by numeri- 
cally solving the Wannier equation (fJJ. 

In the vicinity of the Dirac points |lj], we can write the 
electron-hole dispersion as £k ~ 2hvpk. The attractive 
Coulomb interaction is described by the potential V£ = 
Vk_ p (l ± cos #k.p)/2, calculated with the tight-binding elec- 
tron and hole wave functions in the vicinity of the Dirac 



points, similarly to the approach in Refs. 11281 12911 . Due to 
the finite extension of the p orbitals perpendicular to the 
plane, real graphene is not strictly two-dimensional and the 
2D Coulomb matrix element V^ D has to be modified ac- 
cordingly. In addition, the dielectric environment of graphene 
screens the Coulomb interaction by the effective background 
e. Asa combined effect, one finds 
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Vq = - V™F(qd) , 



(3) 



where the form factor F(qd) depends on the effective 
graphene thickness d. The major effect of F is to supply an 
intrinsic length, needed to fix the length and energy scale. 
As shown in Fig. [T] the Coulomb potential can be approxi- 
mated quite well by an Ohno potential with F(x) = e~ 5x . 
For the numerical evaluations, we take d = 0.18A which is 
obtained from a density functional computation 13611 . Assum- 
ing an s-like ground state, we numerically solve Eq. (0 using 
the linear dispersion and the full graphene Coulomb matrix 
elements. More specifically, we have chosen vf = 9.07 • 10 5 
m/s, yielding «g ~ 2.41 for graphene in vacuum (e = 1) and 
vary the effective background dielectric constant e to change 
«G oc 1/e. 

We find the very large binding energy E\ s = — 3.8 eV 
for the lowest exciton state of graphene in vacuum (e = 1, 
«g = 2.41). This energy rapidly decreases with increasing 
e (decreased a G ), e.g. e = 2 (a c = 1.21), 3 (a c = 0.80), 
and 4 (a Q = 0.60) produce E H = -0.47 eV, -44meV, and 
— 0.94 me V, respectively. As discussed below, these numbers 
scale strictly as 1/rf if the effective thickness of graphene is 
changed. Fig. [2^ shows the full e (lower scale) and aQ (up- 
per scale) dependency of Ei s , Within the numerical accuracy, 
we find a critical value e; on sa 4.8 (ai on « 0.5). For e < £j on 
(ac > a; on ), we always find a negative value for E\ s , indi- 
cating the existence of a bound exciton state. Above (below) 
the threshold value our computations produce strictly positive- 
valued Ei s , indicating that only ionized electron-hole states 
exist. Hence, our results show that the presence of bound 
electron-hole pairs, i.e. strong Coulomb interaction effects, 
depends extremely sensitively on the dielectric environment 
of the graphene. 

The transition from the strongly to a weakly Coulomb in- 
teracting regime can be detected even more clearly by study- 
ing the ratio between the Coulomb and the kinetic energy 
computed using the numerically determined wave functions 
0i s (k). In Fig. |2j5, we plot g = E cou \/Ea n as function of 
aG (lower scale) and e (upper scale). For the bound states, 
the Coulomb energy dominates the kinetic energy so that g is 
larger than one, while the kinetic energy exceeds the Coulomb 
energy for ionized electron-hole pairs. The shaded area in Fig. 
I2J5 indicates the region of ionized states. We see that g crosses 
from the region of bound states to the region of ionized states 
roughly at ej on = 4.8 (ai on = 0.5), in full agreement with 
the ionization threshold observed in Fig.|2^. If we regard g as 
an order parameter, the exciton binding experiences a second- 
order phase transition at a; on because dg/dac becomes dis- 
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FIG. 2. (top) Exciton binding energy vs. background dielectric con- 
stant (lower axis) and qg (upper axis). The left axis presents the 
scaled energies and the right axis gives the energy in eV calcu- 
lated using d = 0.18 A. The full numerical solution (solid line) 
is compared with results obtained using the ansatz wave functions 
f(y) oc (1 + y 2 y 3/2 (dashed line) and f(y) oc e~ M ~ (dotted line). 
The shaded area corresponds to the analytic result l(6}. (bottom) Ratio 
of the Coulomb and kinetic energy vs. ota (lower axis) and e (upper 
axis). The shaded area indicates the region of ionized solutions. 



continuous at ac = ai on - 

In order to analyze the graphene Wannier equation fur- 
ther, we introduce the dimensionless quantities k = kd, 
E\ = E\d/2hvp, 4>(k) = 4>(k)/d, and take the continuum 
limit ^ k — > A/{2tt) 2 Jd 2 k. This yields the dimensionless 
equivalent to Eq. © 



Hx(k) 



4tt 



m*Wl 



k'l) 



|k-k' 



<Mk') = ^A(k) ( 4 ) 



which has to be solved subject to the normalization condi- 
tion J d 2 fc|0A(k)| 2 = 1. Equation (|4]i shows that the effective 
graphene fine-structure constant ac = e 2 /4ireoehvF is the 
single parameter combination characterizing the solution. In 
particular, any solution is independent of the sheet thickness 
d and hence, valid for all graphene-like systems. The effective 
thickness is needed only to fix the length and energy scale. 
Moreover, one recognizes that any solution of Eq. © with 
a finite eigenvalue E\ produces a divergent physical energy 
eigenvalue in the strict 2D limit d — >• 0. 

Assuming rf ^ 0, and using the dimensionless scaling pa- 
rameter u = a/d, we treat the energy eigenvalue E\ as func- 
tional of the scaling parameter: 



E x = 
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with V x {u) = ±Jd 2 yd 2 y'rx(r,u) 
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and TJ™1 = Jd 2 yy n \ijjx(y; u)\ 2 . Stationary points of Eq. © 
correspond to either bound states or resonances of the system, 
depending on the sign of E\ . 

If it exists, the lowest bound state defines the characteristic 
length and energy scales of the exciton, i.e. the exciton Bohr 
radius and binding energy, in terms of the sheet thickness d 
and the energy unit Ihvpjd. Since the wave function scales 
with the sheet thickness which is independent of ac and the 
stationary points of Eq. © depend on ac via the wave equa- 
tion, the scaled wave function ip(y; u) has an implicit nontriv- 
ial u dependence. As a consequence, Tl n ] and the Coulomb 
integral also depend on u. In contrast, the implicit u depen- 
dence of the total energy vanishes for any solution of the wave 
equation. Hence, stationary points of Eq. © are found for 
dE x /du = -Ex/u - acu- 1 dV{u)/du = 0. Since V{u) is 
a monotonically increasing function of u, this predicts stable 
bound states for any value of a with E\ < 0. From these argu- 
ments, we obtain the condition ac > a\ on = T^/Vx(u). For 
every regular, normalized wave function, T'- 1 ' > and V(u) 
does not diverge, and hence, graphene has a finite ioniza- 
tion threshold aj on . This value separates the strong Coulomb- 
interaction regime where bound states exist from the weakly 
interacting configurations without bound states. 

An analytical estimate for aj on can be obtained using a vari- 
ational analysis with a trial wave function ^'(y; it) — > f(y) 
with fixed width to calculate the ratio g = -E C oui/-Ekin of the 
Coulombic and kinetic energy contributions. Figure |2] shows 
the results for f(y) = \/8tt(1 + y 2 )" 3 ^ 2 with the ionization 



threshold aj on = 7r/4 (dashed line) and f(y) = \/87re -y 
with cti on = 1 (dotted line) as function of aG- We notice that 
the resulting ionization values of the fine-structure constant 
depend on the explicit form of the trial wave function. Clearly, 
these estimates only provide an upper bound since the true 
wave function may vary its shape via the u dependence to 
minimize the energy. 

The same analysis can be repeated for a model system 
with a parabolic electron-hole dispersion £t = h 2 'k 2 /2m 
parametrized by the mass m. This situation corresponds to 
a two-band model for an effectively 2D quantum well sys- 
tem. Here, we obtain Eft = -^j [T^ - f-V(d/a)] where 
ciq = 2irh 2 eeo/me 2 assigns a new length scale. Whereas 
this cannot be done for graphene, for the parabolic model 
we may now set d = to proceed to the ideal 2D case. In 
this situation, the kinetic energy scales like 1/a 2 while the 
Coulomb energy is proportional to 1/a. The variational con- 
dition dEft/da = -E/a - E kin /a = -2E/a + E coa] /a = 
yields the well known virial theorem, the Bohr radius a = ciq, 
and the binding energy Eft = — „ h s T^ for any bound state 
solution regardless of e. Hence, there is no ionization thresh- 
old in parabolic -band semiconductors. 

Returning to graphene, we now derive an expression which 
explicitly shows the extreme sensitivity of E\ s on e (or ao). 
For this purpose, we combine the properties resulting from 
Eq. © with the numerical results. A direct differentiation 
of Eq. © produces dE ls /da c = %-^ ~ V{u)/u = 



— Ecoxa/ao since dE\ s /du vanishes for the ground state. As 
suggested by Fig.|2] we can use E coyA /E^ m -l = c(a G -aion)- 
A straightforward integration gives 
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(6) 



where _E™ C is the binding energy of graphene in vacuum and 

P = 1/cCKion. 

Based on Fig. |5] we use P = 4.5, ei on = 4.8262, and 
Elf = —0.0542. In particular, we see that E\ s scales like 
(^ion — e) 4,5 /e — > (eion — e) 4 ' 5 close to the ionization threshold 
which explains the rapid decrease of the binding energy be- 
fore the ionization is reached. The shaded area in Fig.|2]shows 
the analytic estimate of E\ s together with the numerically ob- 
tained binding energy (solid line), nicely demonstrating that 
Eq. © captures the essential features. 

In conclusion, our calculations predict an extreme sensi- 
tivity of the Coulomb effects in graphene on the dielectric 
properties of the environment. We find strongly bound exci- 
tons for freestanding graphene in vacuum, whereas all bound 
states disappear for the effective background dielectric con- 
stant eion = 4.8 which corresponds to the graphene fine- 
structure constant ac = 1/2. Our results are in general agree- 
ment with an analytical real space analysis of the 2D Dirac 
two-body problem 13211 . studies of the Coulomb scattering at 
an impurity charge Ze in g raphene 13 3143511 . as well as Monte 
Carlo calculations Il7lll8ll . predicting instabilities for critical 
a values in the range of 0.5 to 1.66. Physically, this situa- 
tion is realized, e.g., when graphene is deposited onto sub- 
strates or sandwiched between dielectric media such that the 
effective background dielectric constants exceeds 4.8. While 
in the regime of weak Coulomb interaction quasi-free elec- 
trons and holes are responsible for the low-energy electronic 
and optical properties, those are dominated by bound excitons 
for small background dielectric constants. On this basis, we 
expect a pronounced influence of the dielectric environment 
on the graphene ground state and on properties like the low- 
density conductivity, the quantum Hall effect, or the terahertz 
response. 
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